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Abstract 

We present the next-to-leading order QCD corrected rate for the production of a scalar Higgs 
boson with a pair of high pj- bottom and anti-bottom quarks at the Tevatron and at the Large 



Cf^ . Hadron Collider. Results are given for both the Standard Model and the Minimal Supersymmetric 

o _ 

Standard Model. The exclusive bbh production rate is small in the Standard Model, but it can 



be greatly enhanced in the Minimal Supersymmetric Standard Model for large tan /3, making bbh 
^ • an important discovery mode. We find that the next-to-leading order QCD results are much less 

sensitive to the renormalization and factorization scales than the lowest order results, but have a 
^ ' significant dependence on the choice of the renormalization scheme for the bottom quark Yukawa 



coupling. 
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I. INTRODUCTION 



One of the most important problems of particle physics is to uncover the origin of the 
electroweak symmetry breaking. In the simplest version of the Standard Model (SM) of 
particle physics, the breaking of the electroweak symmetry introduces a single physical scalar 
particle, the Higgs boson, that couples to both gauge bosons and fermions. Extensions of 
the Standard Model, like the Minimal Supersymmetric Standard Model (MSSM), introduce 
several scalar and pseudoscalar Higgs bosons. Finding experimental evidence for one or more 
Higgs particles is therefore a major goal of current and future accelerators. Direct searches 
at LEP2 require that the SM Higgs boson mass (M/j) be heavier than 114.4 GeV (at 95% 
c.l.) 1|, while precision electroweak measurements imply < 219 GeV (at 95% c.l.) P|. 
The light scalar Higgs particle of the MSSM (h^) should have mass between the theoretical 
upper bound of about 130 GeV and the experimental lower bound from LEP2, M^o >91 GeV 
(at 95% C.L, 0.5 < tan/3 < 2.4 excluded) j3[. In both Higgs boson should lie in a 

mass region which will certainly be explored at either the Fermilab pp Tevatron collider or 
at the CERN pp Large Hadron Collider (LHC). 

The dominant production mechanism for a SM Higgs boson in hadronic interactions is 
gluon fusion. Among the subleading modes, the associated production with either elec- 
troweak gauge bosons or top quark pairs, as well as weak boson fusion, play crucial roles. 
The inclusion of higher order QCD corrections is in general essential to stabilize the the- 
oretical predictions of the corresponding rates. All of them have now been calculated at 



next-to-leading order (NLO) 



rresponaing rates. All oi tnem nave no\ 

Q Q 1^ iiy, y, y Q y y, y, Q 



and, in the case of 



gluon fusion and associated production with gauge bosons, at next-to-next-to-leading order 
(NNLO) y , llll La y, y, bll in perturbative QCD. 

If the Standard Model is not the full story, however, then other mechanisms of Higgs 
production become very important. Here, we focus on Higgs boson production with a pair 
of bottom quark and antiquark. The coupling of the Higgs boson to a bb pair is suppressed in 
the Standard Model by the small factor, rrib/v, where = (-\/2G'ir)~^/^ = 246 GeV, implying 
that the SM Higgs production rate in association with bottom quarks is very small at both 
the Tevatron and the LHC. In a two Higgs doublet model or in the MSSM, however, this 
coupling grows with the ratio of neutral Higgs boson vacuum expectation values, tan j3, and 
can be significantly enhanced over the Standard Model coupling, leading to an observable 
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production rate for a Higgs boson in association with bottom quarks in some regions of the 
parameter space. 

The production of a Higgs boson in association with bottom quarks at hadron colhders 
has been the subject of much recent theoretical interest. At the tree level, the cross section 
is almost entirely dominated by gg — hbh^ with only a small contribution from qq bbh, 
at both the Tevatron and the LHC The integration over the phase space of the final state 
bottom quarks gives origin to large logarithms proportional to ln(m;,//i/j) (where Hh — Mh), 
which arise from the splitting of an initial gluon into a pair of almost on-shell collinear 
bottom quarks. The use of bottom quark parton distribution functions in the proton (or 
anti-proton) sums these large logarithms to all orders, and could therefore improve a fixed 
order calculation. The inclusive cross section for bbh production should then be dominated by 
the bottom quark fusion process bb —>■ h, as originally proposed in Ref. js^. Some important 
progress has been achieved recent 



NNLO in QCD 23]. At NLO 



y. The bb ^ h production process has been calculated at 
, the residual factorization scale dependence is quite 
large, but at NNLO there is almost no scale dependence. Interestingly enough the NNLO 
results show that the perturbative cross section is better behaved when the factorization 
scale is /ij~M/i/4 (and the renormalization scale is fir — Mh), as expected on quite general 



theoretical grounds 



26 



27 



29|. Moreover, the inclusive bbh cross section has been 



obtained at NLO in QCD via a fixed order calculation that includes the O(a^) corrections 

nnn 

to the parton level processes gg,qq bbh |3fll. 13 IL I32|. The obtained results are compatible 
with hb h dX NNLO, and show that there is actually no large discrepancy between the 
NLO fixed order calculation and the use of 6-quark parton distribution functions, contrary 
to what was originally claimed. However, the results of the fixed order calculation have a 
substantial scale dependence and a better control of the residual large uncertainty is desirable 
for a complete understanding of the comparison between the two approaches. 

In spite of its theoretical interest, the inclusive cross section is experimentally relevant 
only if a Higgs boson can be detected above the background without tagging any of the 
outgoing bottom quarks. Higgs production from bb fusion could be useful, for instance, 
in a supersymmetric model with a lar ge v alue of tan/3, when combined with the decays 



T ' T 



35, Hal. However, even in this case, the 



inclusive measurement of a Higgs signal would not determine the bottom quark Yukawa 
coupling unambiguously, since it should be interpreted as the result of the combined action 
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of other production channels besides bb h^, (e.g. gg — > H^). 

Requiring one or two high pt bottom quarks in the final state reduces the signal cross 
section with respect to bb h, but it also greatly reduces the background Moreover, 
it assures that the detected Higgs boson has been radiated off a bottom or anti-bottom 
quark and the corresponding cross section is therefore unambiguously proportional to the 
bottom quark Yukawa coupling. Using arguments similar to the ones illustrated above 
for the case of the inclusive cross section, one can argue that if the final state has one 
high pt bottom quark then the relevant subprocess is gb ^ bh 24 L The cross section 
for gb — s> bh has been computed including NLO QCD corrections 2S| and the residual 
uncertainty due to higher order QCD corrections is small. On the other hand, if the final 
state has two high pj- bottom quarks and a Higgs boson, then no final state bottom quark 
can originate from a bottom quark parton distribution function. The lowest order relevant 
parton level processes are unambiguously gg bbh and qq — bbh. While the rate for 
this final state is considerably smaller than for the bb ^ h and gb bh subprocesses, the 
background is correspondingly reduced. The final states can be further categorized according 
to the decay of the Higgs boson. Existing studies have considered mostly the dominant 
Higgs decay channel, h —>■ bb 0, [s^, 



h + 



0,Q- 



4G 



41 



44| . but also h t^t [37| and 



In this paper, we present the NLO QCD corrected rates and phase space distributions 
for the fully exclusive processes pp,pp —>■ bbh, where the final state includes two high pt 
bottom quarks. In order to reproduce as closely as possible the currently used experimental 
cuts, we require the final state bottom quark/ ant i- quark to have a transverse momentum 
higher than j9^* = 20 GeV and a pseudorapidity |?7| <2 for the Tevatron and |?7| <2.5 for the 
LHC. The cut on greatly affects the cross section and we therefore study the dependence 
of the cross section on this cut. Similar results have been recently presented in Ref. js^, 
where however no cut on the pseudorapidity has been imposed. Our discussion will focus 
on assessing the uncertainty of the theoretical prediction for the exclusive pp,pp — > bbh 
rates, after the full set of NLO QCD corrections has been included. We will show how the 
large dependence on the unphysical renormalization and factorization scales present in the 
lowest order (LO) calculation of the cross section is greatly reduced at NLO. Moreover, we 
will study the dependence on the choice of renormalization scheme for the bottom quark 
Yukawa coupling. While for Higgs decays and Higgs production in e^e~ collisions using 



4 



the MS definition of tlie bottom quark Yukawa coupling is an efficient way of improving 
tfie perturbative calculation of the corresponding rate by resumming large logarithms at all 
orders 



47 



49 



this may be less compelling in the case of hadronic Higgs production. 
Finally, we will extend our calculation to the scalar sector of the MSSM, including the SM 
QCD corrections at NLO. Preliminary results of the study described in this paper have been 
already presented at several conferences jsil. 

The plan of the paper is as follows. In Section|n]we present an overview of our calculation. 
Since the NLO QCD corrections to gg, qq — » hhh proceed in strict analogy to those for 



qq 



tth 



10 



11 



ij, Q and gg tth |K 



12 



Iq . we will be very brief on details 



and devote more time to the discussion of the residual theoretical uncertainty, emphasizing 
those aspects that are characteristic of the hhh production process. Numerical results for 
the Tevatron and the LHC will be presented in Section 11111 for both the SM Higgs boson 
and the scalar MSSM Higgs bosons in some prototype regions of the model parameter space. 
Section HVl contains our conclusions. 



II. CALCULATION 



A. Basics 



The total cross section for pp,pp hhh at 0{al) can be written as: 
<yNLo{pp,PP bih) = 

^ I - I [j^f{Xi,fi)J^'^{x2, X2, fi) + {l^ 2)] , (1) 

where !Ff'^ are the NLO parton distribution functions (PDFs) for parton i in a proton or 
anti-proton, defined at a generic factorization scale /i/ = /i, and a^^Lo is the 0{a^) parton- 
level total cross section for incoming partons i and j, made of the channels qq, gg hhh and 
{Qj 0)9 ~^ hbh{q, q), and renormalized at an arbitrary scale fir which we also take to be fir = f^- 
Throughout this paper we will always assume the factorization and renormalization scales 
to be equal, fir = fif = fi. The partonic center-of-mass energy squared, s, is given in terms 
of the hadronic center-of-mass energy squared, Sh, by s = xiX2Sh- At both the Tevatron 
and the LHC, the dominant contribution is from the gluon-gluon initial state, although we 
include all initial states. 
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FIG. 1: Sample of diagrams corresponding to 0{as) virtual corrections where the Higgs boson 
couples to an internal fermion loop and not to the external bb pair. The circled cross denotes all 
possible insertion of the final state Higgs boson leg, each insertion corresponding to a different 
diagram. 



The NLO parton-level total cross section reads 



(2) 



where ^lo{xi,X2, fi) is the C(a^) Born cross section, and Scr]^Lo{xi,X2, fi) consists of the 
0{as) corrections to the Born cross sections for gg, qq bbh and of the tree level (g, q)g — >■ 
bbh{q, q) processes, including the effects of mass factorization. 

The evaluation of a^,,n proceeds along the same lines as the corresponding calculation for 



tth production 



11 



12 



13 



14 



15l | and we refer to Refs. 



3 



for a detailed description 



of the techniques used in our calculation. We notice that, in view of the generalization to 
the MSSM with a very enhanced bottom quark Yukawa coupling, both top and bottom 
quark loops need to be considered in those virtual diagrams where the Higgs boson couples 
directly to a closed loop of fermions, a sample of which is illustrated in Fig. ^ 

Contrary to the case of tih production, the NLO cross section for bbh production depends 
significantly on the renormalization scheme used for the bottom quark Yukawa coupling, 
i.e. for the bottom quark mass appearing in g^ij^ = nib/ v. In our calculation of the NLO 
PPi PP ~^ bbh cross section we have considered, for the renormalization of the bottom quark 



Yukawa coupling, both the on-shell and the MS subtraction schemes (in the tth case we only 
used the on-shell top quark renormalized mass everywhere ^^). The MS scheme results in 
a running bottom quark Yukawa coupling and potentially gives better control over higher 
order contributions beyond the 1-loop corrections. We will study the origin and magnitude 
of the residual scheme dependence in Sections IllBI and 1111 Al 
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B. Renormalization scheme dependence 



The ultraviolet (UV) divergences arising from self-energy and vertex 0{as) virtual cor- 
rections to qq, gg —>■ bbh are regularized in = 4 — 2e dimensions and renormalized by in- 
troducing counterterms for the wave functions of the external fields {^^Z^^ (for q = u,d, c, s), 
SZ2^ for the bottom quark, and 6Z^ for the gluon), for the bottom quark mass, 6mb, and 
for the bottom quark Yukawa and strong coupling constants, Sg^i^ and 6Za^ . We follow the 
same renormalization prescription and notation adopted in Refs. jisl . 15 1 for the NLO tih 



inclusive cross section. Consequently, we fix the wave-function renormalization constants 
of the external massless quark fields, SZ2'^\ using on-shell subtraction, while we define the 



wave function renormalization constant of external gluons, SZ^, using the MS subtraction 



scheme and the as renormalization constant, SZ^^, using the MS scheme modified to de- 



couple the top quark 



Refs. 



13 



i3- 



53[. Explicit expressions for 6Z2'^\ 6Z^, and 6Za^ can be found in 



However, given the large sensitivity of the MS bottom quark mass to the renormalization 
scale and given the prominent role it plays in the bbh production cross section through the 
overall bottom quark Yukawa coupling, we investigate here the dependence of the final 
results on the renormalization prescription adopted for the bottom quark. We consider both 



the on-shell (05*) and MS subtraction schemes, for both the bottom quark mass and wave 
function renormalization constants. 

When using the OS subtraction scheme, we fix the wave function renormalization con- 
stant of the external bottom quark field, (5Z2*'')o5, and the mass renormalization constant, 
(5mfc)os, by requiring that 

hip = mfe) = ; lim = , (3) 

where 

S6 = - rUb) (Sv + 5Zf ) + rrib {^s + Sy - (4) 

denotes the renormalized bottom quark self-energy at 1-loop in QCD, expressed in terms of 
the vector, Ey, and scalar, S5, parts of the unrenormalized self-energy, and of the mass and 
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wave function renormalization constants. Using Eq. Q in (i = 4 — 2e dimensions one finds 



OS 47r \ ml J \euv 



where we have exphcitly distinguished between ultraviolet and infrared divergences. The 
infrared divergences are cancelled between virtual and real soft and collinear contributions 



according to the pattern outlined in Refs. 



13,ll5f, to which we refer for more details. 



In the MS scheme, the bottom quark renormalization constants are fixed by requiring 
that they cancel the UV divergent parts of the bottom quark self energy of Eq. i.e. 

^) =_|iC,(4.)<r(l + .)A. (8) 



According to the LSZ prescription j5J], one also needs to consider the insertion of the 
renormalized one-loop self-energy corrections on the external bottom quark legs. While 



these terms are zero in the 05* scheme (see Eq. (jHl)), they are not zero in the MS scheme. 
Together with {6Z2^)jYg, their contribution to the NLO cross section equals the contribution 
of the wave function counterterm in the 05* scheme, {6Z2^^)os, as expected from the LSZ 
prescription itself. The cross section does not depend on the renormalization of the external 
particle wave functions. 

We therefore focus on the scheme dependence induced by the choice of different subtrac- 
tion schemes for the bottom quark mass. We note that the bottom quark mass counterterm 
has to be used twice: once to renormalize the bottom quark mass appearing in internal 
propagators and once to renormalize the bottom quark Yukawa coupling. Indeed, if one 
considers only QCD corrections, the counterterm for the bottom quark Yukawa coupling, 

dgbbh=—^ , (9) 

coincides with the counterterm for the bottom quark mass, since the SM Higgs vacuum 
expectation value v is not renormalized at 1-loop in QCD. This stays true when we generalize 
the gi,bh coupling from the SM to the case of the scalar Higgs bosons of the MSSM. 



At 1-loop order in QCD, the relation between the pole mass, m^, and the MS mass, mb{ij)^ 



is indeed determined by the difference between the OS and MS bottom mass counterterms. 



'^■'':SCT, since 



mb{n) =mb<l 





+ 4 


1 = mfe 









An 

or Eq. 



5CT{fi) 



(10) 



Adopting the 05* or MS prescription consists of using either Eq. (joj) or Eq. (jSj) for the bottom 
mass counterterms while substituting or rnb{^) respectively in both the bottom quark 
propagator and Yukawa coupling. At 0{q^^) the two prescriptions give identical results. 
Indeed, replacing nih by rni,{^) in the Yukawa coupling adds a term 



27r 



-6CT{fi)ar,o + Oiat 



(11) 



to the NLO parton level cross section, which compensates exactly for the difference in the 



OS and MS counterterms. On the other hand, using the MS mass in the bottom quark 
propagator. 



rUb 



a. 



1 + irub — 6CT(fi)- 
Att 1 



nib 



+ 0{al) , 



(12) 



of the LO cross section leads to an extra contribution to the MS NLO cross section which, 
together with the MS mass counterterm insertions into the internal bottom quark propaga- 
tors (see diagrams Si in Fig. 2 of Ref. [l^ and S'2, 6*3, and 5*4 in Fig. 2 of Ref. 13]), coincides 
with the corresponding mass counterterm insertions in the 05* scheme at (9(a^). 



Therefore, using OS or MS at 0{a^g) is perturbatively consistent, the difference between 
the two schemes being of higher order and hence, strictly speaking, part of the theoretical 
uncertainty of the NLO calculation. One notices however that some of the large loga- 
rithms involved in the renormalization procedure of the NLO cross section come from the 



renormalization of the bottom quark mass, and are nicely factored out by using the MS 
bottom mass in the bottom quark Yukawa coupling (see Eq. (fTUj) ). Therefore one should 
consider reorganizing the perturbative expansion in terms of leading logarithms (of the 
form Q;"(/i) ln'^(/i^/m^)) or next-to-leading-logarithms (of the form a"(yu) ln"~^(/i^/m^), for 



/i ~ Mh)-, as obtained by replacing the MS bottom mass in the Yukawa coupling by the 



corresponding 1-loop or 2-loop renormalization group improved MS masses: 



mb{n)ii 
mb{ij)2i 



ITLb 



nib 



OisijUb) 

as{nib) 



1 Co , , asifi) - as{mb) 

l + — {ci- hi) 

60 vr 



1 



4 ag(mfc) \ 

3 TT y 



(13) 



(14) 



where 




) 




1 



(15) 



1 51N - 19nif 
27r UN -2n,f 



(lOliV-lOn^y) , 



(16) 



are the one and two loop coefficients of the QCD /3-function and mass anomalous dimension 
7m, while iV = 3 is the number of colors and n;j = 5 is the number of light flavors. 

In both Higgs boson decays to heavy quarks and Higgs boson production with heavy 
quarks in e'^e~ collisions, using Eq. (fT^ at LO and Eq. (fT^ at NLO in the Yukawa coupling 
proves to be a very powerful way to stabilize the perturbative calculation of the cross sec- 
tion 47]. The difference between LO and NLO rates is reduced and the dependence on the 
renormalization and factorization scales at NLO is very mild, indicating a very small resid- 
ual theoretical error or equivalently a very good convergence of the perturbative expansion 
of the corresponding rate. This is due to the fact that in these cases to a large extent the 
0{as) QCD corrections amount to a renormalization of the heavy quark mass in the Yukawa 
coupling. In more complicated cases, like the case of the hadronic cross section discussed in 
this paper, the previous argument is not automatically true. 

Using the OS or MS bottom quark mass mainly affects the Yukawa coupling. There- 
fore, in the hadronic case, we will look at the different behavior of the NLO cross section 
when the bottom quark Yukawa coupling is renormalized either in the OS or in the MS 
scheme, keeping the bottom quark pole mass everywhere else. Fig. |21 of Section IIIII shows 
the renormalization and factorization scale dependence of the LO and NLO cross sections 
for pp, pp bbh obtained using in the Yukawa coupling either the pole mass or the MS 
running mass rribifJ') Eq. (fT!?|) (at LO) and p4|l (at NLO). The use of mf,(/i) both at LO 
and NLO seems to improve the perturbative calculation of the cross section, since the NLO 
MS cross section is better behaved than the NLO OS cross section at low scales and since 
the difference between LO and NLO cross section is smaller when the bottom quark Yukawa 
coupling is renormalized in the MS scheme than in the 05* scheme. However, both the 
OS and the MS cross sections have very well defined regions of minimum sensitivity to the 
variation of the renormalization/factorization scale and these regions do not quite overlap. 
The difference between the OS and MS results at the plateau should rather be interpreted, 
in the absence of a NNLO calculation, as an upper bound on the theoretical uncertainty. 

The origin of the large difference between the OS and MS NLO cross sections illustrated 
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in Fig. 121 can be understood by studying the numerical effect of the higher order terms that 



bbh {ij = qq, gg) in the OS and MS prescription 



are included in the NLO MS cross section when mf,(/i) is used in the Yukawa coupling. The 
parton level NLO cross sections for ij 
explained above can be written as: 



a 



NLO, OS 



Xl,X2,/i) 



g'io{xi,X2) 



+ 



47r 



mt 



9NLoi^u fi) - 2g'^^{xi,X2)6CT{fi) + —g'J{xi,X2) 



'^JViO.MS"' 



Xi,X2,fi) 



+ 



9Loi^uX2) 



An 



9NLoiXU X2, fJ-) + 



9li{xi,X2) 



(17) 



(18) 



where the dependence on the renormalization scale is explicitly given. Q;s(/i) is the 2-loop 



strong coupling, is the bottom pole mass, and mf,(/i) is the bottom quark MS mass. 
9LOy 9 NLO and gl^ have been defined in such a way that they are the same in the OS 



and the MS schemes. They correspond respectively to the 0{al) {g^io) and 0{al) {g^NLo) 
contributions to the NLO QCD cross section, from which we have singled out the 0{as) 
virtual corrections where the Higgs boson couples to a top quark in a closed fermion loop 
((7*;, see, e.g., diagrams in Fig. H} as well as 6CT{fi), i.e. the difference between the OS and 



MS bottom mass counterterms defined in Eq. (fTUI) . Using Eqs. ()17|) and ()18|) . one can easily 
verify that the difference between the parton level NLO cross sections obtained by using 



either the 05* or the MS scheme for the bottom quark Yukawa coupling is, as expected, of 
higher order in a^, i.e.: 



A 



a 



NLO, OS 
2/ 



0"> 



+ 



(^s{^A9'io{^l,x2) 



m 



2 :^2 



9NLoi^l^^2,f^) + 



nit 



■9li ixi,X2) 



(19) 



The term in the first bracket of Eq. ()19p vanishes at 0{a^), as can be easily verified by using 
Eq. (fTUI) . Hence all the terms in Eq. ()19|) only contribute at 0{a'^) and higher. However, 
while the first term is in general quite small, the term proportional to gNLoi^i, ^2, /^) can be 



large and has a non trivial scale dependence that we can formally write as: 



9 NLO ' 



xi,X2,fi) = gi{xi,X2) + gi{xi,X2)\n [ y 



(20) 
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From renormalization group arguments one can see that g^^{xi,X2) is given by: 



gWxi,X2) = 2 <( (47r&o + 4)^^^o(xi,X2) - / dziPik{zi)g%{xiZi,X2) 

k L-'p 



+ 



dZ2 Pjk{z2)gto{Xl,X2Z2) 



(21) 



where p= {2mh + Mh)'^ I s, Pij{z) denotes the lowest-order regulated Altarelli-Parisi splitting 
function 55] of parton i into parton j, when j carries a fraction z of the momentum of parton 
i, (see e.g. Section V of Ref. and bo is given in Eq. (fT5|) . As a result, A, defined in 

Eq. (fTTHl . turns out to have a non trivial scale dependence and, thus, the difference between 



the NLO hadronic cross section calculated with the OS or with the MS definition of the 
bottom quark Yukawa coupling can be numerically quite significant for some values of the 
renormalization/factorization scale, as we will illustrate in Section UTTl fsee Figs. |2] and El). 



III. NUMERICAL RESULTS 



Our numerical results are obtained using CTEQ5M parton distribution functions for the 
calculation of the NLO cross section, and CTEQ5L parton distribution functions for the 
calculation of the lowest order cross section The NLO (LO) cross section is evaluated 
using the 2-loop (1-loop) evolution of as{fi) with a^^*^(M^) = 0.118. The bottom quark 
pole mass is taken to be mi, = 4.6 GeV. In the OS scheme the bottom quark Yukawa coupling 
is calculated as g^ih = ^b/v, while in the MS scheme as gbbhif^) =^b{f^)/v, where we use 
mi,{fi)u from Eq. (fT^ for a^^o and fnh{fi)2i from Eq. (fT^ for cTjvlo- 

We evaluate the fully exclusive LO and NLO cross sections for bbh production by re- 
quiring that the transverse momentum of both final state bottom and anti-bottom quarks 
be larger than 20 GeV {p^ > 20 GeV), and that their pseudorapidity satisfy the condition 
\r]f,\ <2 for the Tevatron and \rih\ <2.5 for the LHC. This corresponds to an experiment mea- 
suring the Higgs decay products along with two high px bottom quark jets that are clearly 
separated from the beam. Furthermore, we present LO and NLO transverse momentum and 
pseudorapidity distributions. In order to better simulate the detector response, the gluon 
and the bottom/anti-bottom quarks are treated as distinct particles only if the separation 
in the azimuthal angle-pseudorapidity plane is A/2 > 0.4. For smaller values of AR, the four 
momentum vectors of the two particles are combined into an effective bottom/anti-bottom 
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FIG. 2: 0" MLo and a lq for pp bbh at ^/s = 2 TeV (top) and for pp — > at \fs = 14 TeV (bottom) 
as a function of the renormalization/factorization scale /i, for M/^ = 120 GeV. The curves labeled 
Clo.os and (Tjvz,o,os use the OS renormalization scheme for the bottom quark Yukawa coupling, 
while the curves labeled aLo,Ms and (Jnloms use the MS scheme. 
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FIG. 3: The absolute value of the percentage difference A(%) = {(Jmlo.os — (^nlo,'ms)/{'^nlo,os + 
CjvLo.Ms) for pp bbh at ^/s = 2 TeV (top) and for pp — > bbh at ^/s = 14 TeV (bottom) as a 
function of the renormalization/factorization scale /i, for = 120 GeV. The OS and MS labels 
refer to the renormalization scheme chosen for the bottom quark Yukawa coupling. The curves 
labeled as A{mh{fj,)ii) and A(mi,{fj,)2i) use the MS bottom quark Yukawa coupling with the 1-loop 
running mass of Eq. ()13() and the 2-loop running mass of Eq. H14|l. respectively. 
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FIG. 4: (Tmlo.ms and (Ti^o,ms for pp — > hbh at y/s = 2 TeV (top) and for pp — > bhh at ^/s = 14: TeV 
(bottom) as a function of the cut imposed on the final state bottom and anti-bottom transverse 
momentum {p\), for Af/j = 120 GeV and fi = fj.o = mh + Mh/2. 



15 



quark momentum four-vector. 
A. Standard Model Results 

In Fig. 121 we show, for Mh = 120 GeV, the dependence of the LO and NLO cross sec- 
tions for pp bbh at the Tevatron (top) and for pp bbh at the LHC (bottom) on the 
unphysical factorization and renormahzation scale, /i, when using either the OS or the MS 
renormahzation schemes for the bottom quark Yukawa couphng. In both the OS and MS 
schemes the stability of the cross section is greatly improved at NLO, given the much milder 
scale dependence with respect to the corresponding LO cross section. The results presented 
in Fig. 121 are obtained by setting /i = /ir = /i/, i.e. by identifying the renormahzation (/x^) and 
factorization (fif) scales. We have checked that varying them independently does not affect 
the results significantly. By varying the scale fi in the ranges 0.7/io < /i < 4/io (Tevatron) 
and 0.5 fiQ < ji < 8fio (LHC), when using the OS scheme for the bottom quark Yukawa 
coupling, and in the ranges 0.4/io < < 2yUo (Tevatron) and 0.2/io < /i < 2/io (LHC) when 
using the MS scheme, i. e. in the plateau regions, the value of the NLO cross section varies 
by at most 15-20% (where fio = mb + Mh/2) . 

As can be seen in Fig. |2l the cross section calculated with g^^if^ in the MS scheme shows a 
better perturbative behavior, since the difference between a^o and a^Lo is smaller. This is 
in part due to the fact that the LO cross section is calculated using rnb{fJ')u and therefore al- 
ready contains some of the corrections from the renormahzation of the bottom quark Yukawa 
coupling that appear in the NLO cross section as well as at higher order. This observation 
seems to justify the use of mf,(/i)i; at LO and mb{f^)2i at NLO. One also observes that the 
MS NLO cross section is better behaved at low values of the renormalization/factorization 
scales. At the same time, both the OS and MS cross sections show well defined but distinct 
regions of least sensitivity to the renormalization/factorization scale. In both cases this hap- 
pens in the region where the LO and NLO cross section are closer. The variation of the NLO 
cross section with /i about its point of least sensitivity to the renormalization/factorization 
scale is almost the same whether one uses the 05* or MS schemes for the bottom quark 
Yukawa coupling. This indicates that the running of the Yukawa coupling is not the only 
important factor to determine the overall perturbative stability of the NLO cross section. 

As discussed in Section III Bl the numerical difference between the two renormahzation 
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schemes can be significant. This is illustrated in Fig. 01 where we plot the absolute values of 
the relative difference, A = {(yMLo,os — (^nlo,jts)/{<^ NLo,os "I" '^NLo,Ms)i betwccu the hadronic 
cross sections (Jnlcos and (Tmlo^jts at both the Tevatron and the LHC. As discussed in detail 
at the parton level in Section fll Bl (see A defined in Eq. ()19|l ). the difference between the 
two schemes is scale dependent and can be very big for small and large scales. At the LHC, 
the relative difference can be well approximated by A = ^AB with A = ^gNLo/ 9lo and 
B = {1 — (jfib/mbY), where gNLo,Lo correspond to the Qnlo.lo contributions of Eqs. (fT7|l and 
(I18|) calculated at hadron level. For instance, at /i = 0.7/io, A = 0.28 and B = 0.57, while 



at /i = 4/io, A = 0.92 and B = 0.66, which shows that the difference between the MS and 
the OS schemes of the bottom quark is not dominated by the running of the bottom quark 
mass as it would be the case when the majority of the NLO corrections can be absorbed in 
the running of mj,. 

From both the observed similar scale dependence of a^Lo in both schemes and the large 
numerical difference due to the corrections that cannot be absorbed in the running of m^. 



we conclude that the use of the MS bottom quark Yukawa coupling should probably not 
be overemphasized. It is probably a good approximation to take the difference between 
(^NLo.os and ai^LoJTs at their points of least scale sensitivity as an upper bound on the 
theoretical error of the NLO cross section, on top of the uncertainty due to the residual 
scale dependence. This would amount to an additional 15-20% uncertainty arising from the 
dependence on the bottom quark Yukawa coupling renormalization scheme. 

In Fig. m we illustrate the dependence of the exclusive cross section on the px cut imposed 
on the final state bottom and anti-bottom quarks, at both the Tevatron (top) and the LHC 



(bottom). We plot the LO and NLO cross sections obtained using the MS bottom quark 
Yukawa coupling. Reducing the pt cut from 25 GeV to 10 GeV approximately increases 
the cross section by a factor of four. However, as the pt cut is reduced, the theoretical 
calculation of the cross section becomes more unstable, because the integration over the 
phase space of the final state bottom quarks approaches more and more a region of collinear 
singularities. Results without a cut on the transverse momentum of the bottom quarks will 
be presented in a later work j57| (see also Ref. js^)- 

Finally, in Figs. El El and|Hlwe plot the LO and NLO transverse momentum (pj-) and 
pseudorapidity {t]) distributions of the final state particles, the bottom and anti-bottom 
quarks and the Higgs boson, both for the Tevatron and for the LHC. Both LO and NLO 
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FIG. 5: Transverse momentum distributions at LO and NLO of the bottom or anti-bottom quark 
with the largest pT- Shown are the distributions for pp bbh production at ^/s = 2 TeV 

(left) and pp — > bbh production at ^/s = 14 TeV (right) in the SM and using the OS scheme for the 
bottom quark Yukawa coupling. At the Tevatron we choose /x = 2mfe + Mh, while at the LHC we 
choose fj, = 2{2mb + Mh). 

differential cross sections are obtained in the SM and using the OS scheme for the bottom 
quark Yukawa coupling. For the renormalization/factorization scale we choose /i = 2mi,+Mh 
at the Tevatron and /i = 2(2mfe + Mh) at the LHC. These two scales are well within the 
plateau regions where the OS NLO cross sections vary the least with the value of /i. Similar 
results can be obtained using the MS bottom quark Yukawa coupling. 

In Fig. |31we show the LO and NLO pt distributions of the bottom or anti-bottom quark 
with highest px, while Fig. IHl displays the px distributions of the SM Higgs boson. The 
pseudorapidity distributions of the bottom quark and the Higgs boson are shown in Fig. [7| 
and Fig. |H1 respectively. The inclusion of the NLO corrections causes the cross sections to 
be more sharply peaked around low p'^^ and around y]b,h = ^- 

B. MSSM Results 

The rate for hhh production can be significantly enhanced in a supersymmetric model 
with large values of tan (3. In the MSSM, the bottom and top quark couplings to the scalar 



18 




100 120 140 160 180 200 



p!^ [GeV] 



p!^ [GeV] 



FIG. 6: Transverse momentum distributions at LO and NLO of the SM Higgs boson. Shown are 
the distributions for pp bbh production at ^/s = 2 TeV (left) and pp bbh production at 
= 14 TeV (right) in the SM and using the OS scheme for the bottom quark Yukawa couphng. 
At the Tevatron we choose fi = 2rrn, + Mh, while at the LHC we choose /U = 2(2m-5 + Mh). 
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FIG. 7: Pseudorapidity distributions at LO and NLO of the bottom quark. Shown are the % 
distributions for pp bbh production at ^/s = 2 TeV (left) and pp bbh production at ^/s = 14 TeV 
(right) in the SM and using the OS scheme for the bottom quark Yukawa coupling. At the Tevatron 
we choose fj, = 2mb + Mh, while at the LHC we choose fi = 2{2mb + Mh). 
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FIG. 8: Pseudorapidity distributions at LO and NLO of the SM Higgs boson. Shown are the r]fi 
distributions for pp bbh production at ^/s = 2 TeV (left) and pp bbh production at ^/s = 14 TeV 
(right) in the SM and using the OS scheme for the bottom quark Yukawa couphng. At the Tevatron 
we choose fi = 2mi, + M^, while at the LHC we choose fi = 2(2mj, + M^). 



Higgs bosons are given by: 



sma - „ cos a 

bbh'' : -g,,ih UN' : ^—^Quh 



cos (3 sin /5 ' 

cos a -^^n sin a 

hbH^ : -g,j,^ ttH^ : —-Quh 

cos p sm p 



where gi,ifi and guh are the SM bottom and top quark Yukawa couphngs, and are the 
hghter and heavier neutral scalars of the MSSM, and a is the angle which diagonalizes the 
neutral scalar Higgs mass matrix 5^. By replacing the SM top and bottom quark Yukawa 
couplings with the corresponding MSSM ones, our calculation can then be straightforwardly 
generalized to the case of the scalar Higgs bosons of the MSSM. The bottom quark Yukawa 
coupling to the MSSM pseudoscalar Higgs boson, A^, is also enhanced at large tan/3. The 
corresponding cross section for bbA^ production can be obtained from our calculation in the 
mi, limit, which we do not consider in this paper. We will present, however, complete 
results for bbA^ production, i. e. for non-zero mi,, in a future study. 



The MSSM Higgs boson masses and t 



le mixing angle a have been computed up to two- 



loop order using the program FeynHiggs 5^. In Tables HI and ITTl we provide the values of the 
input parameters ((M/jO, tan/?) or (Mj^o, tan/5)) and the resulting values of a used in the 
calculation of the top and bottom quark Yukawa couplings to the light and heavy neutral 
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FIG. 9: 

(^NLo.Ms for pp — > bbh production at y/s — 2 TeV (top) and pp — > bbh production at 
^/s = U TeV (bottom) in the SM and in the MSSM with tan/? = 10, 20, and 40. For the Tevatron 
we considered pp bbh^ with Mf^o = 100, 110, 120, and 130 GeV, while for the LHC we considered 
pp blH^ with M^^o = 120,200,400,600, and 800 GeV. For each (M;,o,tan/?) and (Mji^o, tan/?) 
point, the corresponding values of a and Ma are hsted in Tables HI and ITll 
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tan P = 10 


Mho [GeV] 


100 


110 


120 


130 


Ma [GeV] 


102.42 


113.86 


127.95 


264.72 


a [rad] 


-1.3249 


-1.1963 


-0.9054 


-0.1463 


tan /? = 20 


Mfto [GeV] 


100 


110 


120 


130 


Ma [GeV] 


100.61 


110.95 


121.89 


146.72 


a [rad] 


-1.4420 


-1.3707 


-1.1856 


-0.3108 


tan /? = 40 


M;,o [GeV] 


100 


110 


120 


130 


Ma [GeV] 


100.15 


110.23 


120.46 


133.71 


a [rad] 


-1.5007 


-1.4601 


-1.3444 


-0.4999 



TABLE I: Values of a and Ma, computed up to two-loop order by using the program FeynHiggs j59|, 
corresponding to different choices of tan /3 and Mf^o . In the calculation of a and Ma we choose the 
genuine SUSY input parameters as follows: Mg = M^^ = M^^ = M^^ = M^^ = 1 TeV, M^^^ = 2 TeV, 
Ai, = At = M/-^ + cot p, and ^ = M2 = 200 GeV. 

MSSM scalar Higgs bosons. This choice of MSSM parameters takes into account present 
experimental limits on the MSSM parameter space, but represents otherwise just one among 
many possible realizations of the MSSM parameter space. The results obtained with this 
choice of MSSM input parameters illustrate the typical enhancements over the SM results 
one can expect when considering the production of neutral scalar Higgs bosons in association 
with bottom quarks. 

The top part of Fig. El compares the NLO pp — > bbh SM cross section at the Tevatron 
with the corresponding cross section for production of the lightest neutral scalar Higgs 
boson in the MSSM for tan/5 = 10, 20, and 40. A large enhancement of up to three orders 
of magnitude is observed. As the light neutral Higgs boson mass approaches its maximum 
value, the mixing angle a becomes very small, as can be clearly seen in Table HI This has 
the effect of suppressing the bbh^ rates at this point. A similar effect can be observed in 
the production of a heavy neutral Higgs boson when Mho is approaching its minimum value 
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tan (3 = 


10 






Mho [GeV] 


120 


200 


400 


600 


800 


Ma [GeV] 


108.05 


198.55 


399.41 


599.64 


799.74 


a [rad] 


-0.9018 


-0.1762 


-0.1140 


-0.1057 


-0.1030 






tan P = 


20 






Mho [GeV] 


120 


200 


400 


600 


800 


Ma [GeV] 


116.45 


199.56 


399.81 


599.89 


799.91 


a [rad] 


-0.5785 


-0.0901 


-0.0574 


-0.0531 


-0.0517 






tan P = 


40 






Mj^o [GeV] 


120 


200 


400 


600 


800 


Ma [GeV] 


118.92 


199.82 


399.92 


599.95 


799.96 


a [rad] 


-0.3116 


-0.0460 


-0.0289 


-0.0267 


-0.0259 



II: Values of a and Ma, computed up to two-loop order by using the program Feyn- 
2, corresponding to different choices of tan f3 and Mho- In the calculation of a and Ma we 



TABLI 

Higgs 

choose the genuine SUSY input parameters as follows: Mg = M^^ = M^^ = Mj^^ = Mj^^ = 1 TeV, 
M[^ = 0, Ab = At = M/^^ + /icot (3, and ^ = Ms = 1 TeV. 



(see Table ITT]) . as shown in the bottom part of Fig. El Again, we compare the production of 
the SM Higgs boson with that of the heavier neutral scalar Higgs boson of the MSSM and 
observe a significant enhancement of the rate in the MSSM for large tan f3. 



IV. CONCLUSIONS 

We presented results for the next-to-leading order QCD cross section for exclusive bbh 
production at both the Tevatron and the LHC. Our NLO results show an improved sta- 
bility with respect to the unphysical factorization and renormalization scales as compared 
to the leading order results and increase the reliability of the theoretical prediction. The 
uncertainty in the resummation of large logarithms from higher order corrections, however, 
is also visible in the dependence of the NLO cross section on the renormalization scheme of 
the bottom quark Yukawa coupling. The residual renormalization/factorization scale depen- 



23 



dence is of the order of 15-20% when the bottom quark Yukawa couphng is renormahzed in 
the OS or MS schemes respectively. We conservatively estimate the additional uncertainty 
due to the renormalization scheme dependence of the bottom quark Yukawa coupling to be 
at most of order 15-20%. 

Our calculation is important for Higgs boson searches at hadron colliders where two high 
Pt bottom quarks are tagged in the final state. In supersymmetric models with large tan f3, 
bbh production can be an important discovery channel, at both the Tevatron and the LHC 
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